Measurement of absolute transition frequencies of 87 Rb to nS and nD Rydberg states 
by means of electromagnetically induced transparency 
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We report the measurement of absolute excitation frequencies of 87 Rb to nS and nD Rydberg 
states. The Rydberg transition frequencies are obtained by observing electromagnetically induced 
transparency on a rubidium vapor cell. The accuracy of the measurement of each state is < 1 MHz, 
which is achieved by frequency stabilizing the two diode lasers employed for the spectroscopy to 
a frequency comb and a frequency comb calibrated wavelength meter, respectively. Based on the 
spectroscopic data we determine the quantum defects of 87 Rb, and compare it with previous mea- 
surements on 85 Rb. We determine the ionization frequency from the 5Si/2{F—l) ground state of 
87 Rb to 1010.0291646(3) THz, providing the binding energy of the ground state with an accuracy 
improved by two orders of magnitude. 
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I. INTRODUCTION 

Rydberg atoms are attracting large research interest as 
their quantum state can be prepared and controlled with 
high precision and flexibility by means of electromagnetic 
fields Actual research subjects include cavity quan- 
tum electrodynamics [2|, Rydberg- Rydberg interactions 
3] with possible applications in quantum information 
processing [1] , ultracold chemistry Q @ Q > and sensing 
dispersion forces between atoms and surfaces @ [9| [Toj . 
Due to their large polarizability, Rydberg atoms respond 
to static and dynamic electric fields with level shifts, and 
are therefore sensitive to interactions with microscopic 
and macroscopic objects. For precision measurements of 
such interactions and accurate quantum state control, it 
is beneficial to know the unperturbed Rydberg energy 
levels that we quantify in the present paper for the 87 Rb 
atom, which is widely used in experiments. 

In alkali atoms, the deviation from the hydrogen en- 
ergy structure can be expressed in terms of quantum de- 
fects (5 [l[ . If the quantum defects are known, the energy 
of Rydberg levels can be calculated with respect to the 
ionization limit. The quantum defects of the nS and nD 
lines of 85 Rb have been measured by Li et al. us- 
ing microwave excitation between Rydberg states. This 
method provides superior frequency resolution on the or- 
der of ~ 10 kHz, however, it does not serve as an absolute 
frequency reference. If, in addition, absolute transition 
frequencies between the ground and Rydberg states are 
known, the binding energy of the ground state can be 
determined. 

Absolute measurements of the nS and nD lines of 
85 Rb as well as the isotope shift for 87 Rb have been done 
by Stoicheff and Weinberger [l2| using two-photon spec- 
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troscopy with an accuracy of 100 MHz. Higher accuracy 
measurements include the nF states of the 85 Rb isotope 
with an accuracy of 8 MHz [l3| and on the nS levels of 
85 Rb with an uncertainty greater than 6 MHz [l4j]. In the 
present work, we measure Rydberg transition frequencies 
of 87 Rb with an absolute accuracy of < 1 MHz, which re- 
sults from frequency stabilizing the lasers employed for 
the spectroscopy to absolute frequency references. 

High resolution spectroscopy of Rydberg states is pos- 
sible by obser ving electromagnetically induced trans- 
parency (EIT) [15|. The detection of Rydberg states by 
means of EIT, and the measurement of the fine struc- 
ture splitting of 85 Rb have been reported by Mohapatra 
et al. 16]. In the present work, we implement EIT for 
measuring the absolute frequency of Rydberg excitations 
of 87 Rb, from which we determine the quantum defects 
of the nS and nD levels and the ground state ionization 
energy of this isotope. 



II. EIT LEVEL SCHEME AND 
EXPERIMENTAL SETUP 

The energy level diagram of 87 Rb relevant for this work 
is shown in Fig.[TJ The three level system consists of the 
87 Rb 5S 1/2 (F=2) ground state, the 5P 3/2 (F=3) inter- 
mediate state, and nS or nD Rydberg states with a high 
principal quantum number n. We observe EIT spectra 
on a rubidium vapor cell using two counterpropagating 
laser beams. The frequency of the "probe laser" is fixed 
at the 5S 1 / 2 (F=2) <->■ 5P 3 / 2 (F=3) transition, which is 
known to be 384.2281152 THz Q3 (= 780.2460209 nm), 
while the "coupling laser" (ss 480 nm) is scanned across 
the Rydberg resonance. We record EIT spectra by mon- 
itoring the probe laser transmission as a function of the 
coupling laser frequency. 

The experimental setup is illustrated in Fig. [5] Abso- 
lute frequency measurements are achieved by referencing 
the probe and coupling laser frequencies to a frequency 
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FIG. 1. Energy level diagram of Rb and corresponding 
laser wavelengths for observation of electromagnetically in- 
duced transparency using a Rydberg state. The probe laser 
is stabilized to the lower atomic resonance, while the cou- 
pling laser, with a variable detuning hA c , is scanned across 
the Rydberg resonance. 
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FIG. 2. Optical setup (schematic) for the observation of EIT 
in a rubidium vapor cell. The probe and coupling beams 
are counterpropagating in the rubidium vapor cell (heated 
to ~ 45 °C) and the probe laser intensity is monitored with 
a photodiode (PD). A chopper in the coupling laser beam 
provides a reference signal for lock-in detection of the pho- 
todiode signal. The probe laser is frequency stabilized to a 
frequency comb. The coupling beam at 480 nm comes from a 
frequency doubled 960 nm laser that is stabilized and scanned 
by a wavelength meter. The wavelength meter is calibrated 
with the help of a frequency comb beat signal at 960 nm. 



comb (Menlo Systems "FC 1500"). 

The source of the probe beam is a grating stabilized 
diode laser of ~ 500 kHz linewidth. The laser is stabilized 
to 780.2460209 nm by superimposing it with the beam of 
the frequency comb on a fast photodiode and locking 
the corresponding radio frequency beat signal with the 
nearest comb mode. The beat signal deviation, and thus 
the accuracy of the probe laser frequency, was maintained 
within 500 kHz. 

The coupling beam at 480 nm is sourced from a fre- 
quency doubled, grating stabilized diode laser (Toptica 
"TA SHG Pro" ) at 960 nm of - 500 kHz linewidth. In or- 
der to access any desired coupling frequency, we stabilize 
the 960 nm laser to a Fizeau interferometer based wave- 
length meter (HighFinesse "WS Ultimate-2" ) . The wave- 
length meter with a built-in servo feedback unit measures 
frequency, and controls the frequency scans of the cou- 
pling laser. 

In order to assure the absolute accuracy of the wave- 
length meter, we perform an automated calibration pro- 
cedure. A beat signal between the 960 nm laser and a 
known frequency comb mode is recorded with a digital 
spectrum analyzer. The beat frequency is determined 
by software peak recognition and the corresponding laser 
frequency is used to calibrate the wavelength meter. The 
resulting accuracy of frequency control by the calibrated 
wavelength meter was characterized with the frequency 
comb by exemplary measurements of the frequency offset 
with respect to the set frequency. We found the frequency 
offset after calibration approximately to be normally dis- 
tributed with a standard deviation of 400 kHz. 

The experimental sequence for recording EIT spectra 
on Rydberg levels is as follows. First, the wavelength me- 
ter is calibrated using the method described above. Sec- 
ond, a trigger is sent to the wavelength meter, that starts 
sweeping the coupling laser frequency (at most 12 MHz/s 
sweeping rate, 30 MHz span at ~ 960 nm), while record- 
ing the measured coupling laser frequency. The photo- 



diode signal of the probe beam is recorded by a digital 
storage oscilloscope. If the scan range covers an EIT res- 
onance, a sharp transmission peak is recorded. In order 
to detect weak EIT signals of high n states, the coupling 
laser beam is periodically blocked by a chopper wheel at 
a frequency of about 1 kHz, and the photodiode signal 
of the probe beam is recorded using a lock-in amplifier 
synchronized to the chopper. As both the laser frequency 
and the photodiode signal are recorded simultaneously, 
we obtain the probe laser transmission as a function of 
the coupling laser frequency as shown in Fig. [3] The mea- 
surement at each EIT resonance is fitted by a Gaussian 
function to determine the peak position. In order to re- 
duce the statistical error, the measurement sequence is 
performed five times for each Rydberg level. 




(4.(2ttTHz) 

FIG. 3. (Color online) EIT resonance for the 34S , 1/2 Ryd- 
berg state (exemplary) as recorded by the photodiode, against 
the coupling laser frequency measured simultaneously by the 
wavelength meter. A Gaussian shape is fitted for peak fre- 
quency determination (solid line). 
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III. RYDBERG TRANSITION FREQUENCIES 
AND MEASUREMENT ACCURACY 

Absolute frequencies of transitions 5P3/2(-F=3) — > 
nS 1/2 (F=2), 5P 3/2 (F=3) -> nD 3/2 and 5P 3/2 (F=3) -> 
nD 5 / 2 of 87 Rb were observed for principal quantum num- 
bers n in the range 19 — 65 (nS states) and 19 — 57 (nD 
states). The results are listed in tables Hill IIVI andlVl (Ap- 
pendix). We note that we were able to observe Rydberg 
states up to n = 180 with the help of the lock-in am- 
plifier. Such highly excited Rydberg states, however, ex- 
hibit large linewidths and partially irregular lineshapes, 
thus these are not analyzed in the present manuscript. 
For quantum defect calculations, we use only the data 
given in tables Hill HV1 and IV1 

The uncertainty of the frequency data comprise sev- 
eral known error sources of physical and technical na- 
ture. The main error source arises from the frequency 
stabilization of the probe laser with an uncertainty of 
500 kHz (see section HI]) . A frequency deviation A p 
from the atomic resonance lj p leads to the excitation 
of atoms with a velocity of v = c— . If A p is small 
compared to the Doppler broadening, the EIT transmis- 
sion peak appears at a coupling laser frequency shifted 
by A c = ^uj c — (uj c /ujp)Ap with respect to the reso- 
nance frequency u> c . Consequently, the uncertainty of 
the probe laser frequency translates into an uncertainty 
of (780 nm/480 nm) ■ 500 kHz w 800 kHz of the EIT peak 
position measured with the coupling laser. 

The uncertainty of the coupling laser frequency enters 
the EIT peak determination with 200 kHz. This is less 
than the uncertainty of the frequency stabilization of this 
laser (400 kHz) because of the five- fold repetition of the 
frequency measurement (including calibration). 

In our setup, the frequency of the comb modes can 
be determined with an absolute accuracy on the or- 
der of 10 kHz. This uncertainty is negligible with re- 
spect to other sources of error. The high accuracy of 
the frequency comb is achieved by a 10 MHz GPS disci- 
plined rubidium frequency reference (Precision Test Sys- 
tems "GPS10RBN") with a specified Allan deviation of 
2 x 10~ n at Is observation time. The accuracy of the 
transfer of the rf reference frequency to optical frequen- 
cies was determined by Kubina et al [l8[ , who compared 
two similar frequency combs referenced to one common rf 
frequency standard, finding the optical frequencies that 
were measured with both combs to agree to a level of 
6 x 10~ 16 . 

The width of the observed EIT resonances is typi- 
cally < 10 MHz. The linewidth of the lasers, technical 
noise, triggering precision, and physical deviations from 
the Gaussian line shape that has been used for finding 
the peak position, contribute to the uncertainty of the 
measured resonance frequency as well. The overall un- 
certainty resulting from these effects was determined in- 
dependently for each transition, and was typically less 
than 300 kHz. 



Despite the strong n 7 dependence of the polarizability 
[l|, line shifts due to stray electric fields are not present 
in the vapor cell due to screening b y io ns and electrons 
on the inner surface of the glass cell [16| . We verified the 
screening in our setup for field strengths up to 470 V/cm 
and on Rydberg states up to n = 150. We observed line 
broadening due to external magnetic fields, which are in 
our setup on the order of the earth magnetic field, but 
no significant line shift or asymmetry of the resonance 
peak. Pressure shifts at a cell temperature of < 45 °C 
are calculated to be less than 200 kHz [l9[ . 

IV. QUANTUM DEFECTS AND IONIZATION 
ENERGY 

The energies of Rydberg levels [l[ are given by 



where Ei is the ionization energy threshold and 
TV = ^^Koo = h ■ 3289.82119466(2) THz 

+ m S7 Rb 

is the Rydberg constant, corrected for the reduced elec- 
tron mass in 87 Rb. 

For sufficiently large principal quantum numbers n, the 
quantum defects 8(n, I, j) depend only little on n and can 
be approximated by the modified Rydbcrg-Ritz parame- 
ters 

8(n,l,j)nS + - h (2) 
(n - d y 

Within the accuracy of our measurements on n > 19 
Rydberg states, the inclusion of higher order terms in 
the approximation did not lead to improved results. 

A precise measurement of the modified Rydberg-Ritz 
parameters for the nS, nP and nD series of 85 Rb was 
performed by Li et al. [ll| . By exciting microwave tran- 
sitions between Rydberg levels (n — 32 up to 37) and 
subsequent detection by field ionization, the transition 
frequencies were determined with a resolution of up to 
10 kHz. The "average" quantum defects 5(n,n + 1) de- 
termined by this measurement were used to improve the 
value of quantum defects previously obtained by Loren- 
zen and Niemax (2pj . 

Traditionally, the quantum defects are determined for 
the fine structure levels. The hyperfine splitting is neg- 
ligible for nD states, but not for nS states because of a 
nonzero probability amplitude at the position of the nu- 
cleus. In our measurements of nS states, only transitions 
to F=2 states were observed because of selection rules. 
In order to obtain quantum defects for the fine structure, 
we need to consider the hyperfine shifts in the evaluation 
of our data. 

The hyperfine splitting of states between n = 28 and 
33 of both 85 Rb and 87 Rb was determined by Li et al [TT[ | . 
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Using their data, we extrapolate the hyperfine splitting 
between the nSi/ 2 (F=l) and (F=2) states of 87 Rb ac- 
cording to a n*~ 3 law (n* := n — 5(n, l,j)) as 

Ahfs,f= 2 - A H fs,f=i = 33.5(9) GHz • n*~ 3 . 

In general, the hyperfine shift is given by Ahfs,f = 
|(F(F + 1) -7(7 + 1)- J(J + 1)), which, compared with 
the measurement, yields the value of the hyperfine con- 
stant A. For the J = 1/2 states of 87 Rb (nuclear spin / = 
3/2), the hyperfine splitting is Arts f=2 ~ Ahfs f=i — 
2A. 

Extending the extrapolation with the n*~ 3 law for the 
states n > 28 to n = 19 leads to a systematical error 
in the determination of A. In order to quantify this er- 
ror we extrapolate further down to n = 11, for which 
the hyperfine splitting is known from a measurement 
by Farley et al. [2l|. We find 10% deviation between 
the calculated (A±is 1/2 = 34.5 MHz) and the measured 
(Ans 1/2 = 37.4(3) MHz) value, that gives an upper limit 
of the systematic error for n > 11 states. 

We now can compute the hyperfine shift of the 
nS\/ 2 (F=2) states that have been measured in this work, 

A H fs,f=2 = \ A =Y7l- 34 ( 3 ) GHz ' n *~*> 

and subtract these frequency shifts from the measured 
transition frequencies to the nS\/ 2 (F=2) states to find 
the quantum defects with respect to the fine structure 
levels. For the lowest nS state measured in this work, 
n = 19, the hyperfine shift is Ahfs,f=2 = 3.2(3) MHz. 

In order to determine the modified Rydberg-Ritz pa- 
rameters, equation (UJ) is simultaneously fitted to all three 
data series (nSx/2, nD 3 / 2 an d nD 5 / 2 ), with the ioniza- 
tion frequency uii = E\/% as a common fit parameter. 
The difference between measured frequencies and corre- 
sponding values obtained by the fit (residuals) are shown 
in Fig. 2] exemplary for the measurement series on the 
nSi/2 lines. The residuals vary around zero with an RMS 
deviation of < 0.75 MHz, which applies to all three data 
series, consistent with the expected overall precision of 
the experimental setup. 

The resulting values for the modified Rydberg-Ritz pa- 
rameters are shown in table U along the values determined 
by Li et al. (ill for 85 Rb. The parameters for the nS se- 
ries of 85 Rb [ll| and 87 Rb (this work) match each other 
with comparable uncertainties. Omitting the hyperfine 
structure correction, the fit of our data returns a slightly 
different value of So = 3.1311788(8) (with no difference 
in #2), that gives the value of the quantum defect of the 
nSi/ 2 (F=2) hyperfine levels. 

The values of So of the nD series show a deviation from 
the values of Li et al. of about three times the stated 
uncertainty. Nonetheless, comparing the frequency dif- 
ferences measured in this work (u c ,n+i — w Ci „) with the 
Rydberg - Rydberg transition frequencies w n , n +i mea- 
sured by Li et al., we find that the values available in 
both data sets are well within one standard deviation 
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FIG. 4. Difference between measured transition frequencies 
to nSi/2 Rydberg states (see table [TTT) ) and the transition fre- 
quencies calculated with equation ([T} based on the Rydberg- 
Ritz parameters and ionization frequency specified in tables |TJ 
andUJl) 

(nSi/2 and nD 5 / 2 ) or at most just outside one standard 
deviation (nZ? 3 / 2 ). We therefore suspect differences in 
the analysis as well as the range of included n values to 
be the reason for the discrepancy. 

The value of the ground state ionization energy (i.e. 
the binding energy of the ground state), as obtained by 
fitting equation (TTJ to the experimental data, is given in 
table mi Because all measured Rydberg transition fre- 
quencies enter into this value, its statistical error is small 
(90 kHz) and systematic uncertainties like pressure shifts 
become significant. We estimate systematic shifts to be 
on the order of 200 kHz, resulting in < 300 kHz uncer- 
tainty of the ground state ionization energy. The ground 
state ionization energy found in this work is within the 
uncertainty stated by Stoicheff and Weinberger [l^ and 
is two orders of magnitude more accurate. 







85 Rb (Li et al [11]) 


87 Rb (this work) 


nS 1/2 


So 


3.131 180 4(10) 


3.131 180 7(8) 




82 


0.178 4(6) 


0.178 7(2) 


nD 3/2 


So 


1.348 091 7(4) 


1.348 094 8(11) 




5 2 


-0.602 9(3) 


-0.605 4(4) 


nD 5/2 


So 


1.346 465 7(3) 


1.346 462 2(11) 




82 


-0.596 0(2) 


-0.594 0(4) 



TABLE I. Modified Rydberg-Ritz parameters 80 and 82 of 
the nS and nD lines of 85 Rb following and of 87 Rb as de- 
termined in this work. nS quantum defects refer to the fine 
structure levels, corrected for shifts due to hyperfine struc- 
ture. 



V. CONCLUSIONS 

In summary, we performed precision spectroscopy of 
Rydberg states of 87 Rb using diode lasers, frequency sta- 
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Stoicheff and 
Weinberger [12J 


this work 


Ei,5P 3/2 (F = 3)/h 




625.794 214 8(3) 
THz 


£i,5S 1/2 (F=l)/(7!.c) 
Ei,5S 1/2 (F=l)/h 


33 690.945 5(15) 

cm -1 
1 010.029 14(5) 

THz 


33 690.946 44(1) 

cm" 1 
1 010.029 164 6(3) 
THz 



TABLE II. Ionization frequency from the 5Ps/2{F=3) state 
(fit result including systematic error) and the 5Si/2(F=l) 
ground state (calculated on basis of the fit result) of 87 Rb as 
determined in this work and compared with the value given 
by Stoicheff and Weinberger [l2| . 



bilized to a frequency comb and a calibrated wavelength 
meter. This allows for absolute frequency measurements 
with an uncertainty of < 1 MHz. The system presented 
in this work combines the ease of use of the wavelength 
meter with the accuracy of the frequency comb, enabling 
absolute frequency measurements over a wide frequency 
range, not limited by the availability of atomic lines to 
serve as an absolute reference. 

We determined quantum defects and the ionization en- 
ergy of 87 Rb with an uncertainty of < 300 kHz. In the 
future, the precision of such measurements can be im- 
proved further by reducing the laser linewidths and in- 
creasing the bandwidth of the frequency locks. 
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Appendix: Rydberg transition frequencies 
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TABLE III: Measured absolute transition frequencies (lj c ) 
between the 5Ps/2(F=3) state and nSi/2(F=2) Rydberg 
states of 87 Rb. The quantum defects 8 include the correc- 
tion for hyperfine structure. Aco c is the difference between 
the measured transition frequency and the frequency cal- 
culated with eq. JTJ using the quantum defects and ion- 
ization energy determined in this work. 

2. nD 3/2 



n 


(2tt THz) 


6 


Acj c 
(2tt MHz) 


19 


615 


2383571 


8 




1 


3461520( 




-0.1 


20 


616 


3395780 


'8 




1 


3463543( 


l\ 


0.3 


21 


617 


2770769 


8 




1 


3465269( 


10 




0.3 


22 


618 


0817641 






1 


3466765 ( 


11 




-0.9 


23 


618 


7775900 


» 




1 


3468033( 


13 




0.1 


24 


619 


3833361 


'8 




1 


3469144(14' 




0.3 


25 


619 


9139034 


8 




1 


347013 


1 


) 


-0.4 


26 


620 


3812369 


'8 




1 


347100 






-0.4 


27 


620 


7949991 


'8 




1 


347174 


2 




0.2 


28 


621 


1630743 


8 




1 


347241 


"'2 




0.5 


29 


621 


4919489 


'8 




1 


347301 


'3 




0.6 


30 


621 


7869975 


'8 




1 


347356 


3 




0.3 


31 


622 


0527045 


8 




1 


347407 


3 




-0.2 


32 


622 


2928357 


'8 




1 


347452 


'4 




-0.3 


33 


622 


5105718 


'8 




1 


347488 


V 4 




0.6 


34 


622 


7086101 


8 




1 


347533 


} 4 




-1.1 


35 


622 


8892621 


'8 




1 


347556 


'5 




0.8 


37 


623 


2060305 






1 


347618 






0.1 


38 


623 


3453297 


; 




1 


347640 




[ 


0.6 


39 


623 


4736776 


8 




1 


347668 


7 




-0.0 


40 


623 


5921942 


'8 




1 


347689 


'7 




0.1 


41 


623 


7018578 


8 




1 


347711 


'8 




-0.1 



6 



42 
43 
44 
45 
46 
47 
48 
51 
53 
54 
55 
57 
TABLE 



623 
623 
623 
624 
624 
624 
624 
624 
624 
624 
624 
624 

Tvr 

5P 3/2 (F=3) 



(2tt THz) 



8035293(8) 
8979667(8 
9858378(8 
0677426(8 
1442056(8 
2156986(8 
2826438(8 
4597844(8 
5611212(8 
6075157(8 
6513410(8 
7320082(8 
Measured 
-¥ n_D 3/2 of 



347727(8) 

347740(9) 

347778(10) 

347778(10 

347789(11 

347808(12) 

347827(13) 

34784(2) 

34788(2) 

34788(2 

34786(2 

34789(2 



Acj c 
(2tt MHz) 



0.1 
0.5 
-1.4 
-0.1 
0.2 
-0.3 
-0.7 
0.7 
-0.7 
-0.2 
0.9 
0.4 



iydberg transition frequencies 
7 Rb. See description of table 
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5/2 



n 


(2tt THz) 
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(2tt MHz) 


19 


615.2402660 


8) 
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0.6 


20 
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8) 


1.3447550 
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-0.1 


21 


617.2784664 


8) 


1.3449236 


9) 
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618.0829639 


8) 


1.3450698 
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-0.2 
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8) 
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619.3842463 


81 


1.3453060 


15 
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619.9147046(8) 
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345402(2) 


-0.7 


26 


620.3819464(8 
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27 
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345560(2) 


-0.2 


28 


621.1636344(8) 
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345629(2) 
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29 
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30 
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1 


345738(3) 
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1.0 
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1 


345901(4) 
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37 
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1 
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0.5 


38 
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0.1 


40 
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1 
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0.9 


41 
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1 
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42 
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-0.1 


43 
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1 


346112(9) 


0.7 


44 
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1 
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346162(10 


-0.9 


46 
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1 
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1 
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51 
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34622(2) 
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53 


624.5611995(8 


1 


34624(2) 


-0.2 


54 


624.6075890(8 


1 


34626(2) 


-0.4 


55 


624.6514104(8) 


1 


34623(2) 


0.9 


57 

T71 1 


624.7320704(8) 


1 


34626(2) 


0.5 



Measured Rydbcrg transition frequencies 
5^3/2(^=3) — > nD 5 / 2 of 87 Rb. See description of table 
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